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I  INTRODUCTION 


The  following  discussion  highlights  the  need  to  characterize 
simultaneously  the  electrochemical  and  mechanical  states  in  a  growing 
stress  corrosion  or  corrosion-fatigue  crack.  We  also  point  out  the  need 
for  a  suitable  method  for  investigating  electrochemical  behavior  within 
a  crack.  The  research  program  described  here  addresses  these  needs  and 
ultimately  aims  at  using  the  resulting  data  to  develop  a  predictive 
model  for  stress  corrosion  cracking  (SCC)  and  corrosion  fatigue  (CF). 

The  load-bearing  capacity  of  most  alloys  can  be  weakened  by  a 
corrosive  environment.  Two  common  forms  of  environment-sensitive 
mechanical  property  degradation  are  SCC  and  CF.  Both  forms  greatly 
influence  the  safety,  reliability,  and  economics  of  many  components  used 
in  military  and  commercial  power  generation,  transportation,  and  energy 
conversion  systems.  Therefore,  substantial  effort  has  been  made  to 
understand  environment-assisted  crack  growth,  to  identify  and  develop 
materials  that  have  Improved  environment  resistance,  and  to  formulate 
computational,  predictive  models. *-*^ 

As  a  result,  various  models  have  been  proposed  over  the  years  for 
SCC,  CF,  or  other  types  of  environment-assisted  crack  growth.  Wei  and 
Landes,*®  for  example,  considered  the  crack  growth  rate  in  a  CF 
situation  to  consist  of  the  sum  of  the  crack  growth  rate  in  an  inert 
environment  and  an  environment-assisted  growth  rate.  The  Wei  and  Landes 
theory  has  been  shown  to  explain  adequately  many  features  of  the  CF 
behavior  of  Ti  and  T1-6A1-AV  in  NaCl  solutions.**  Furthermore,  this 
work  emphasized  the  need  for  a  better  understanding  of  the  processes 
that  occur  at  the  crack  tip,  particularly  at  threshold  load  levels  for 
crack  growth  and  at  low  loading  frequencies  (<  1  Hz). 

For  SCC  alone,  Logan*'  proposed  a  model  for  cracking  in  aqueous 
environments  based  on  a  slip-dissolution  mechanism  at  the  crack  tip. 

This  mechanism,  now  widely  accepted, relates  crack  growth  to  the 
enhanced  anodic  dissolution  current  at  the  crack  tip  resulting  from  the 
intermittent  mechanical  rupture  of  an  environmentally  created  film  at 
the  crack  tip. 
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Logan's  model  has  been  modified  to  include  a  variety  of  rate¬ 
controlling  processes  that  include  repassivaton  rate,  microcreep, 

o 

solution  renewal,  and  critical  charge  concepts.  If  the  processes 
occurring  within  the  crack  were  better  understood,  it  would  be  possible 
to  determine  the  appropriate  slip-dissolution  model  for  any  given 
material,  environment,  and  load  application. 

Quantitative  modeling  of  the  cracking  process  is  complicated 
because  it  Involves  complex  interactions  between  the  cracking  material, 
the  environment,  and  the  applied  load.  Because  specimen  microstructure, 
surface  chemistry,  electrochemistry,  and  loading  conditions  are  expected 
to  Influence  cracking  behavior,  they  should  be  accounted  for  in  attempts 
to  model  crack  growth.  Clearly,  simultaneous  characterization  of  the 
mechanical  and  electrochemical  state  at  the  crack  tip  is  required  for 
development  of  a  reliable,  predictive  crack  growth  model. 

The  mechanical  state  at  the  crack  tip  is  usually  described  in  these 
models  using  classical  fracture  mechanics  parameters  such  as  the  linear 
elastic  stress  intensity  (K  parameter)1^*17  or,  more  recently,  the 

IQ  IQ 

elastic-plastic  stress  intensity  (J  parameter),  *  .  Unfortunately, 

the  electrochemical  state  of  the  crack  tip  is  not  so  easily  character¬ 
ized  because,  until  recently,  no  reliable  experimental  method  existed 
that  could  either  provide  measurements  directly  within  the  crack  or 
could  convert  measurements  made  external  to  the  crack  into  crack  tip 
parameters.  Some  attempts  have  been  made  to  measure  the  local 
chemistries  that  develop  at  crack  tips  by  using  microanalytlcal  chemical 
techniques  on  specimens  that  are  actually  cracking  or  on  artlflcal 
crevices  that  had  been  constructed  to  simulate  stress  corrosion 
cracks.  However,  most  of  past  crack  growth  models  have  used 
electrochemical  data  obtained  on  bulk  specimens  and  have  assumed  that 
such  data,  such  as  the  free  corrosion  potential,  are  representative  of 
values  inside  the  crack,  or  these  bulk  values  were  extrapolated  to  crack 
tip  values  by  means  of  some  postulated  model.  As  might  be  expected, 

crack  growth  rates  predicted  by  such  models  are  several  orders  of 

20  21 

magnitude  different  from  those  that  are  actually  observed.  * 


One  possible  method  for  evaluating  the  electrochemical  state  of 
material  in  the  vicinity  of  a  crack  tip  is  the  AC  impedance  technique 
This  technique,  under  development  for  the  past  several  years,  ha: 

only  recently  been  applied  to  the  problem  of  environmentally  assisted 
cracking.  This  approach  shows  promise  for  determining  crack  tip 
electrochemical  parameters  by  applying  the  transmission  line  model  to 
deconvolve  the  measured  AC  impedance  spectra  of  growing  cracks.  A 
detailed  description  of  this  technique  as  applied  to  the  problem  of 
environmentally  assisted  cracking  is  given  in  the  following  section. 


II  THEORETICAL  CONSIDERATIONS 


The  stress  corrosion  or  corrosion-fatigue  crack  Is  modeled  by  a 
cracked  electrode.  The  solid  and  electrolyte  phases  in  the  cracked 
electrode  are  represented  by  resistance  elements  parallel  to  the 
direction  of  the  corrosion  current  flux,  and  the  interfacial  resistance 
is  represented  by  an  impedance  distributed  between  the  other  two 
resistances.  The  electrode  may  be  treated  as  ''one-dimensional,,  if  the 
potential  (E),  current  flux  (I),  and  reactant  concentrations  (c)  vary 
only  with  the  depth  within  the  crack  and  not  with  the  lateral  position 
along  the  crack.  In  this  case,  the  local  values  of  E,  I,  and  c  may  be 
replaced  by  their  average  values  in  a  line  perpendicular  to  the  crack 
plane. 
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Ksenzheck  and  Stender  discussed  the  conditions  that  are 
applicable  to  transforming  a  three-dimensional  problem  into  a  one¬ 
dimensional  one.  The  mathematical  model  for  the  electrical 
representation  shown  in  Figure  1  requires  solutions  to  a  set  of 
differential  equations  of  the  form 
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where  r  is  the  resistance;  subscripts  e  and  s  refer  to  the  electrolyte 
and  solid  phases,  I  ■  I  +  Ig,  Z  is  the  lnterfaclal  impedance,  x  is 
distance  into  the  crack,  and  AE  ■  E  -  E  . 


Electrolyte 


FIGURE  1  TRANSMISSION  LINE  MODEL  FOR  A  CRACK  IN  AN  AQUEOUS  ENVIRONMENT 


For  the  simplest  set  of  boundary  conditions  corresponding  to  an 
electrolyte-filled  right  cylindrical  crack  in  which  rg  and  re  are 
independent  of  x  and  the  solid  and  electrolyte  phase  current  contacts 
are  diametrically  opposed,  as  shown  in  Figure  1,  the  transmission  line 


impedance  Zt^  is  given  by 
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where  1  is  the  current  crack  length,  and 

Y  -  f(rs  +  re)/Z1]1/2  (5) 

The  interfacial  impedance,  Zj,  in  equation  (5)  is  equivalent  to  that  at 
a  plane-parallel  electrode.  That  is,  it  contains  the  Information  that 
would  be  available  if  it  were  possible  to  determine  the  electrochemical 
kinetic  properties  of  each  infinitesimal  element,  dx,  of  the  electrode 
thickness  in  the  x  direction. 

22—24 

The  conventional  equivalent  circuit  for  Z^  is  shown  in  Figure  2 

The  components  of  this  circuit  contain  the  electrochemical  kinetic 
parameters  that  are  important  in  evaluating  the  behavior  of  a  cracked 
electrode.  Brief  descriptions  of  these  parameters  follow. 

Double  Layer  Capacitance,  C.  The  double-layer  capacitance,  C,  is 
directly  proportional  to  the  area  of  the  electrolyte/solid  phase 
Interface.  Thus,  given  a  knowledge  of  the  double-layer  capacitance  per 
unit  area,  C',  under  the  prevailing  conditions  of  reactant  concentra¬ 
tion,  temperature,  and  potential,  one  can  determine  C  by  deconvolving 

the  measured  transmission  line  Impedance.  This  value  of  C  can  then  be 

28 

used  to  estimate  the  total  wetted  area  per  unit  electrode  thickness, 
which  we  define  as  A'.  Values  of  C'  may  be  available  from  the 
literature  or  can  be  measured  directly  for  a  plane-parallel  electrode  of 
known  area. 


FIGURE  2 


Double  Layer 
Capacitance 

JA-350583-29 

INTERFACIAL  IMPEDANCE  ELEMENT; 
RANDLES'  EQUIVALENT  CIRCUIT 


Warburg  Coefficient,  a.  At  moderate  frequencies  (>  1  Hz),  the 
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Warburg  impedance,  Zw,  can  be  expressed  as 


Zw  -  a(l  -  j)w1/2 

where 


(6) 


o 


RT 

n2F2cA"(2D)1/2 


(7) 


c  is  the  reactant  concentration,  and  A"  expresses  the  electrode 
electrolyte  area  per  unit  length  that  is  actively  engaged  in  reaction. 
The  other  parameters  in  equation  (7)  have  their  usual  meanings.  The 
relationship  between  A'  and  A”  can  be  used  to  define  an  area  efficiency 
that  accounts  for  the  presence  of  passive  zones,  which  do  not  contribute 
to  the  cell  current  during  crack  growth. 

Diffusion  Layer  Thickness,  6.  At  the  limit  of  low  frequency,  an 

apparent  resistive  shunt  to  the  Warburg  impedance  is  observed. 

This  resistance  defines  the  steady-state  DC  operating  behavior,  and  its 

value  may  be  used  in  conjunction  with  a  to  calculate  the  effective 

diffusion  layer  thickness,  6.  In  the  derivation  of  equation  (6),  semi- 

infinite  linear  diffusion  was  assumed.  However,  in  the  wide  regions 

near  the  crack  mouth,  natural  convection  will  limit  the  thickness  of  the 

diffusion  layer,  whereas  in  the  narrower  regions  near  the  crack  tip,  the 

dimension  of  the  diffusion  layer  is  proportional  to  the  crack  opening. 

22  30-35 

Thus,  we  can  account  for  the  steady-state  DC  behavior  ’  in  the 

model  of  the  cracked  electrode  by  introducing  the  Nernst  diffusion  layer 
thickness  as  a  boundary  condition. 

Interfacial  Reaction  Resistance,  Rf.  The  reaction  resistance 

together  with  A”  and  6  can  be  used  to  obtain  Information  regarding  the 
14  29  36 

electrode  kinetics  *  *  within  the  cracked  electrode.  The  exchange 

current  density,  iQ,  can  be  determined  from  the  relationship 


R_  “  RT/nFA'i_ 


(8) 


When  used  in  conjunction,  with  the  Stern-Geary  relationship,  '  Rf 
provides  a  measure  of  the  instantaneous  corrosion  rate  within  the  crack. 

Electrolyte  Resistance,  rfi.  The  electrolyte  resistance  per  unit 
length  in  the  cracked  electrode,  re,  gives  a  measure  of  the  amount, 
composition,  and  distribution  of  electrolyte  in  the  crack.  The 
electrolyte  resistance  will  be  controlled  primarily  by  the  crack  width 
and  to  a  smaller  extent  by  the  diffusion  of  reactant  species  from  the 
bulk  electrolyte  and  the  diffusion  of  reaction  products  from  the  crack 
tip. 

Solid  Electrode  Resistance,  rfi.  The  solid  phase  resistance  per 
unit  length  in  the  cracked  electrode,  rg,  provides  a  measure  of  the 
resistive  component  in  the  cracked  electrode. 

The  reduction  of  the  Impedance  data  to  obtain  the  parameters  of 
interest  for  cracked  electrodes  Involves  applying  nonlinear  regression 
analysis  to  an  Impedance  model  or  equivalent  circuit.  Broers  and 
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coworkersJO  have  demonstrated  the  applicability  of  the  transmission 
line  model  for  characterizing  porous  electrodes  in  molten  carbonate  fuel 
cells,  and  we  believe  that  these  methods  can  be  applied  to  a  growing  SC 
or  CF  crack. 

An  additional  impedance  element,  Zt^p,  has  been  Included  in  Figure  1 
to  account  for  the  special  electrochemical  characteristics  of  the  crack 
tip.  For  example,  if  the  crack  grows  by  a  mechanism  of  film  rupture  and 
anodic  dissolution  at  the  crack  tip,  then  Ztjp  can  be  represented  as  a 
special  case  of  the  equivalent  circuit  in  Figure  2,  with  R,.  determined 
by  the  exchange  current  density  for  metal  dissolution  (possibly 
involving  terms  for  precipitation  of  oxide  and  repassivation).  Under 
the  conditions  of  a  cyclic  mechanical  load,  the  area  term  A' ,  in 
equation  (8)  must  also  be  given  as  a  periodic  function  about  a  mean 

9 

value,  Aq.  In  the  simplest  case,  A'  would  be  linear  function  of  the 
mechanical  load  for  crack-opening  displacement.  Thus, 


where  w  is  the  frequency  of  the  load  variation,  and  $  is  the  phase 
delay  between  load  variation  and  film  rupture. 


Equation  (9)  suggests  a  second  method  of  evaluating  the 
electrochemistry  of  the  crack  tip  environment,  using  the  transmission 
line  model.  Under  a  steady  electrical  bias  (a  small  DC  voltage  more 
positive  than  the  free  corrosion  potential),  a  DC  current  can  be  made  to 
flow  from  the  specimen  to  a  suitable  counter  electrode.  This  current 
originates  with  the  anodic  dissolution  at  the  crack  tip  and  is  given  by: 

I  (u>t)  -  r)/Rr  (10) 

where  T]  ■  V  appHe(|  -  V  free  corrosion* 

Substituting  Equation  (9)  into  Equation  (8),  and  then  substituting 
the  resulting  expression  for  Rr  into  Equation  (10)  yields: 
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RT 
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(10' 


This  expression  for  I(wt)  describes  the  crack  tip  current  expected  for 
cyclic  load  conditions. 

This  can  be  viewed  as  a  sinusoidal  current  source  originating  at 
the  crack  tip.  The  current  source  is  situated  where  the  interfacial 
reaction  resistance  (Rj.  in  Figure  2)  is  placed  in  Zt^p  in  Figure  1. 

Alternatively,  we  can  define  an  electromechanical  impedance, 

zem  *  Z/1  <U> 

which  may  be  examined  as  a  function  of  the  frequency  (and  amplitude)  of 
the  mechanical  load  variation. 

The  advantage  of  an  electromechanical  impedance  study  is  that  the 
perturbation  is  applied  exactly  at  the  point  of  Interest  (the  crack 
tip).  The  alternative  method,  a  variable  frequency  electrical  potential 


variation  applied  to  the  mouth  of  a  crack  in  a  specimen  cycled  with 
constant  frequency  mechanical  load,  results  in  a  perturbation  that 
decreases  in  amplitude  at  the  crack  tip,  as  the  crack  Increases  in 
length. 


Ill  EXPERIMENTAL  METHODS 


Background 

To  calculate  accurately  the  many  electrochemical  parameters  that 
can  be  obtained  from  AC  measurements,  one  must  make  high  precision 
Impedance  measurements  over  a  wide  range  of  frequencies,  particularly  at 
very  low  frequencies.  The  unavailability  of  adequate  Instrumentation 
has  been  the  principal  limitation  to  the  characterization  of  cracked  or 
porous  systems  by  AC  techniques  and  Is  largely  responsible  for  the  only 
modest  success  of  Broers  and  coworkers. ^8-40  ^ore  recently,  Mund  et 

ii 

al.  analyzed  the  performance  of  hydrogen  electrodes  In  alkaline  fuel 
cells  systems  from  AC  Impedance  measurements  applied  to  a  transmission 
line  model. 

The  current  work  of  Hills  and  coworkers^-^  provides  good 
verification  for  the  application  of  the  transmission  line  model  and  the 
extension  of  equilibrium  AC  measurements  for  predicting  the  steady-state 
electrochemical  parameters.  In  one  set  of  experiments,  the  interfacial 
impedance  components  were  calculated  from  the  measured  Impedance  of  a 
packed  bed  of  graphite  spheres  using  the  transmission  line  model,  and 
the  results  were  compared  with  the  interfaclal  Impedance  measured 
directly  for  a  single  sphere. 

Figure  3  shows  the  impedance  data  obtained  for  a  packed-bed  elec¬ 
trode  superimposed  on  the  least-squares  best  fit  of  data  to  equation  (4) 
The  tabulated  values  of  the  double-layer  capacitance,  Warburg  coef¬ 
ficient,  and  exchange  current  density  for  the  packed  bed  have  been 
calculated  from  this  regression  fit.  These  values  may  be  compared  with 
the  Interfaclal  impedance  parameters  measured  directly  for  a  single 
graphite  sphere  Isolated  from  the  bed.  The  substantial  agreement 
between  the  two  sets  of  calculated  impedance  parameters  is  strong 
verification  for  the  applicability  of  the  transmission  line  model  to 
characterizing  porous  or  cracked  electrodes. 
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FIGURE  3 

MEASURED  AND  CALCULATED  IMPEDANCE  SPECTRA 
FOR  A  BED  OF  600-xm  GRAPHITE  SPHERES 

It  is  clear  from  the  foregoing  discussion  that  many  important 
operating  parameters  of  porous  and  cracked  electrodes  are  contained 
within  the  measured  impedance  dispersion.  Two  factors  are  fundamental 
in  the  application  of  this  method  to  the  in  situ  determination  of 
electrode  characteristics.  First,  Impedance  measurements  must  be  of 
sufficient  precision  (better  than  0.1%  for  each  impedance  component) 
and  second,  the  measurements  must  be  made  over  a  frequency  range  wide 
enough  (10-^  to  10^  Hz)  to  adequately  resolve  the  impedance  spectrum  and 
to  allow  subsequent  deconvolution  to  obtain  the  electrochemical 
parameters. 

The  requirements  for  a  high  speed,  high  precision,  Impedance 
measurement  system  capable  of  operating  at  very  low  frequencies  have 
been  presented  in  detail  by  Macdonald  and  McKubre.  In  general,  these 
requirements  are  in  conflict.  Maximum  acquisition  speed  is  available 
from  noise  or  transform  methods,  but  these  have  low  precision. ^  High 
precision  is  best  achieved  with  a  frequency-by-frequency  measurement 
technique  such  as  phase-sensitive  detection,  but  the  measurements  are 
time  consuming.  Fortunately,  this  is  not  a  serious  hindrance  when  such 
equipment  is  automated  and  the  experiments  are  run  under  computer 
control. 

Experimental  Methods  for  Current  Program 

Given  the  requirements  discussed  above,  we  developed  an 
experimental  facility  designed  to  obtain  data  on  crack  electrochemistry 
during  CF  and  SCC  of  various  materials  in  aqueous  solutions  at  room 
temperature  and  pressure,  Figure  4  shows  schematically  the  experiment 
arranged  for  Impedance  measurement  of  specimens  under  mechanical  load 
conditions.  Photographs  of  the  setup  are  shown  in  Figure  5 
(electronics).  Figure  6  (test  chamber),  and  Figure  7  (test  specimen). 

As  detailed  in  Figure  4,  impedance  measurements  were  made  using  a 
Solartron  Model  1250  frequency  response  analyzer  (FRA)  operated  under 
microcomputer  control.  The  Model  1250  is  an  automated,  digitally 
demodulated,  stepped  frequency,  impedance  meter  capable  of  0.01% 


precision  for  each  impedance  component  in  the  frequency  range  10  to 
6  x  10*  Hz.  This  Instrument  and  the  preceding  series  of  frequency 
response  analyzers  (Solartron  Model  1172  and  1174)  have  been  used 
extensively  in  recent  years  to  measure  Impedances  in  electrochemical 
systems. ^5-49 

The  FRA  was  operated  under  the  control  of  an  Apple  II-Plus  micro¬ 
computer  with  64  K-bytes  of  Internal  random  access  memory  (RAM)  and  320 
K-bytes  of  external  RAM  (Axlon  "Ramdisk  320”),  equipped  with  a  20  M-byte 
rigid  disk  (Corvus  model  20  MB/M).  The  control  program  was  written  at 
SRI  in  Basic  and  6502  Machine  languages,  and  is  included  as  Appendix  A. 

Some  experiments  were  performed  under  potentiostatlc  control,  using 
a  Princeton  Applied  Research  Model  173/276  potentlostat;  however,  most 
of  the  experiments  were  performed  under  pseudo-gal vanostatlc  control. 

In  the  latter  mode,  a  resistance  several  orders  of  magnitude  larger  than 
the  cell  impedance  is  placed  in  series  with  the  counter  electrode,  and  a 
large  AC  voltage  is  applied.  To  a  good  approximation,  the  current  is 
determined  by  the  AC  voltage  divided  by  the  series  resistance. 
Independent  of  cell  impedance. 

Because  of  difficulties  in  interfacing  our  data  acquisition  system 
with  the  servohydraulic  mechanical  test  facility,  no  data  have  yet  been 
obtained  on  the  electromechanical  impedance  (as  described  in  Section  II) 
The  data  reported  in  Section  IV  are  for  the  results  of  sinusoidal 
current  and  voltage  perturbations  under  static  and  cycled  load 
conditions. 


For  experiments  in  which  the  load  is  cycled,  it  is  necessary  to 
synchronize  the  periods  of  the  mechanical  and  electrical  cycles.  This 
was  achieved  by  integrating  the  Impedance  data  over  one  complete  cycle 
of  the  mechanical  load.  Thus  the  results  represent  an  average  over  a 
full  cycle  of  crack-opening  displacement.  If  these  cycles  are  not 
synchronized,  then  the  data  may  contain  a  large  amount  of  systematic 
error  and  will  inevitably  shew  considerable  scatter.  For  example,  if 
the  impedance  measurements  are  integrated  over  t.  time  that  is  short  with 
respect  to  the  period  of  the  mechanical  load  cycle,  then  multiple 


determinations  (at  different  electrical  perturbation  frequencies)  may  be 
made  during  a  single  mechanical  cycle;  these  will  cover  the  range  from 
maximum  to  minimum  crack  opening.  If  the  two  periods  are  not 
synchonized,  this  will  appear  as  scatter  in  the  data  set. 

Because  it  is  necessary  to  integrate  over  at  least  one  cycle  at  the 
electrical  perturbation,  some  error  will  be  introduced  at  low  frequen¬ 
cies  when  the  Impedance  data  are  integrated  over  more  than  one  (but  not 
an  integer  number  of)  mechanical  cycles.  This  is  inevitable  and  results 
in  some  scatter  in  our  data  at  very  low  frequencies  (<  0.01  Hz). 

An  environmental  chamber  was  constructed  that  attaches  to  a  servo- 
hydraulic  mechanical  testing  machine.  This  chamber  and  the  associated 
attachments  allow  us  to  circulate  aerated  or  deaerated  3.5%  NaCl 
solution  (or  other  electrolytes)  at  controlled  temperatures  past  our 
fatigue  precracked  compact  tension  specimens.  The  crack  growth  can  be 
measured  by  optical  methods  or  by  the  compliance  method,  by  means  of  a 
water-immersible  LVDT  mounted  on  the  specimen.  A  special  composite 
reference/counter  electrode  incorporates  a  silver/silver  chloride 
electrode  in  a  0.25-inch-diameter  Lugin  probe,  positioned  about  0.1  inch 
from  the  start  of  the  fatigue  precrack.  The  fatigue  precrack  is  at  the 
bottom  of  a  2.00-inch-deep  notch  (1.50  inches  wide  by  0.010  inch  thick), 
and  the  tip  of  the  0.25-inch-diameter  Lugin  probe  is  centered  at  the 
bottom  of  this  notch  to  minimize  "edge  effects”  on  the  electrochemical 
measurements.  The  counter  electrode  is  located  in  a  chamber  to  which 
the  Lugin  probe  is  attached.  This  arrangement  allows  us  to  polarize 
only  the  area  at  the  Lugin  probe  tip  and  thus  measure  the  response  of 
this  area. 
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IV  PROGRESS  DURING  THE  FIRST  YEAR 


Introduction 

In  the  current  program,  we  are  applying  the  AC  impedance  technique 
to  determine  the  distribution  of  -’lectrochemical  processes  in  a  static 
and  growing  stress  corrosion  or  corrosion-fatigue  crack.  The  program 
consists  of  three  research  tasks: 

(1)  Build  up  the  facilities  necessary  to  measure  the  AC  impedance 
spectra  of  growing  and  static  stress  corrosion  and  corrosion- 
fatigue  cracks. 

(2)  Measure  the  AC  impedance  spectra  of  several  alloys  in  aqueous 
environments  as  a  function  of  loading,  crack  length,  and 
potential.  Measure  the  AC  impedance  spectra  on  strained  or 
scratched  alloys  if  necessary. 

(3)  Begin  to  deconvolve  the  AC  impedance  spectra  to  electrochemical 
parameters.  Use  the  electrochemical  parameters  to  model  the 
environment-assisted  crack  growth  mechanism. 

During  the  first  year  of  this  program,  we  have  completed  the 
literature  survey,  selected  the  material/environment  systems,  designed 
and  constructed  the  crack  growth  test  facility,  and  performed  several 
initial  experiments.  The  facility  development  was  described  in  Section 
III.  This  section  describes  the  results  of  our  preliminary  experiments. 

We  selected  HY80  steel  and  T1-6A1-4V  alloys  in  a  3.5%  NaCl  solution 
to  be  investigated  using  the  AC  impedance  technique.  These  two 
materials  were  selected  for  the  following  reasons. 

•  Both  materials  have  had  their  SCC  and  CF  resistance  documented 
for  chloride-containing  environments.  Neither  of  these  alloys 
will  easily  undergo  SCC  in  chloride-containing  environments, 
but  both  exhibit  CF  behavior.  However,  the  environment  can 
easily  be  manipulated  to  bring  about  SCC  in  both  alloys. 
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•  The  HY80  and  T1-6A1-4V  should  have  very  different  passivation 
rates,  surface  film  properties,  and  active  (bare  metal) 
corrosion  properties  in  chloride-containing  environments. 

These  differences  in  properties  should  be  reflected  in  the  AC 
impedance  spectra  and  make  their  interpretation  easier. 

•  Both  materials  have  been  used  or  have  been  considered  for  use 
in  marine  applications.  Data  are  available  that  describe  their 
mechanical,  electrochemical,  SCC  and  CF  cracking  behavior  in 
both  synthetic  and  natural  seawater. 

Results 

Figures  8  and  9  show  impedance  data  for  the  growth  of  a  crack  in 
T1-6A1-4V.  The  data  designated  as  "no  crack"  in  Figure  8  show  the 
impedance  response  in  the  range  10~^  <  f  <  10^  Hz  for  a  specimen 
subjected  to  a  5  Hz  cyclic  load  but  showing  no  superficial  indication  of 
having  cracked.  In  the  Hyquist  domain  this  data  set  appears  as  a 
semicircle  whose  center  is  located  above  the  real  axis.  The  low 
frequency  Intercept  on  the  real  axis  is  very  close  to  2000  Q.  The  high 
frequency  region  is  shown  in  a  considerably  amplified  form  in  Figure 
8(b).  A  second  semicircular  region  is  observed  with  a  high  frequency 
intercept  on  the  real  axis  at  52.8  Q.  The  high  frequency  intercept 
corresponds  to  the  series  resistance  due  to  the  electrolyte  path  between 
the  tip  of  the  reference  electrode  and  the  base  of  the  notch  in  the 
specimen. 

In  the  Bode  plot  [Figure  8(c)],  these  data  exhibit  a  single  minimum 
in  the  phase  occurring  at  approximately  0.15  Hz.  The  magnitude 
information  shows  two  limiting  plateaus  (corresponding  to  the  intercepts 
with  the  real  axis  in  the  Nyquist  plot)  separated  by  linear  regions  of 
slope  -1/2,  -1,  and  -1/4  with  increasing  frequency. 

Examining  these  linear  regions  in  order  of  decreasing  frequency, 
one  sees  that  the  region  at  slope  -1/4  corresponds  to  the  high  frequency 
semicircle  in  the  Nyquist  plot,  shown  in  Figure  8(b).  This  region  is 
present  for  unstressed  specimens  before  cyclic  loading,  and  appears, 
completely  unchanged,  during  the  entire  growth  of  the  crack.  We 
speculate  that  this  region  is  associated  with  diffusion  in  the  plane  of 
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FIGURE  8  IMPEDANCE  RESPONSE  OF  THE  Ti-6A8-4V  SPECIMEN  BEFORE 
AND  AFTER  THE  APPEARANCE  OF  A  MACROSCOPIC  CRACK 
IN  AERATED  3.5%  NaC8  SOLUTION 
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FIGURE  8  IMPEDANCE  RESPONSE  OF  THE  Ti-6A8-4V  SPECIMEN  BEFORE 
AND  AFTER  THE  APPEARANCE  OF  A  MACROSCOPIC  CRACK 
IN  AERATED  3.5%  NaCS  SOLUTION  (Continued) 
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FIGURE  9  THE  EFFECT  OF  LOAD  ON  THE  IMPEDANCE  RESPONSE 
OF  A  MATURE  CRACK  IN  Ti-64E-4V  IN  AERATED 
3.5%  NaCE  SOLUTION 


the  machined  notch  and  perpendicular  to  the  axis  of  the  reference 
electrode.  Although  the  region  of  slope  -1/4  has  significant  extent  in 
the  Bode  plot,  the  Impedances  Involved  are  very  small  and  have  very 
little  significance  in  a  linear  (e.g.,  Nyquist)  plot  of  the  data.  We 
will  not  consider  this  zone  further  in  this  report. 

At  intermediate  frequencies,  all  data  sets  obtained  show  a  region 
of  slope  -1 .  This  corresponds  to  a  f-*  dependence  and  is  the  form 
expected  for  the  equivalent  circuit  in  Figure  2,  at  intermediate  and 
high  frequencies,  when  diffusion  of  metal  ions  from  the  electrode 
surface  are  not  rate  limiting.  A  region  of  slope  -1  in  a  Bode  plot 
corresponds  to  a  semicircle  with  its  center  located  on  the  real  axis  in 
the  Nyquist  representation.  The  beginnings  of  a  semicircle  of  great 
magnitude  (on  the  scale  shown)  can  be  seen  in  Figure  8(b);  the 
parameters  of  this  semicircle  are  determined  by  the  lnterfaclal  reaction 
resistance,  Ry,  and  the  double-layer  capacitance,  C  (see  Section  II). 

Impedance  data  obtained  from  unstressed  specimens  before  the 
application  of  a  cyclic  load  display  a  f~*  dependence  (Bode  slope  *  -1) 
down  to  limltingly  low  frequencies.  Only  under  cyclic  loading 
conditions  does  a  region  of  slope  -1/2  appear,  and  this  region  becomes 
increasingly  significant  as  a  crack  is  observed  to  grow.  Thus,  we 
associate  the  region  of  slope  -1/2  in  a  Bode  plot  with  the  presence  of  a 
crack. 

Figures  8(a),  (b)  and  (c)  each  show  two  curves:  one  labeled  "crack” 
and  the  other  labeled  "no  crack"  on  the  basis  of  visual  inspection  of 
the  sides  of  the  compact  tension  specimen.  Because  both  show  a  region 
of  f~^Z  dependence  in  Figure  8(c),  it  is  apparent  that  the  curve 
designated  as  "no  crack"  has  at  least  a  vestigial  rupture,  associated 
with  cyclic  loading  of  the  oxide  film  at  the  crack  tip.  The  curve 
labeled  "crack”  was  obtained  when  a  crack  of  about  3  mm  length  was 
observed  at  the  side  of  the  specimen. 

“1  /2 

A  region  of  apparent  f  '  dependence  (slope  -1/2)  for  the 
impedance  within  the  crack  is  predicted  by  equations  (4)  and  (5).  The 

Y 1 

transmission  line  impedance,  Ztl»  contains  terms  in  e'  ,  which  can  be 
expanded  as  a  Taylor  series. 


1  +  yl  +  (yl)  /2!  +  (yl)  /3!  +  .  . 


(12) 


For  yl  <  <  1,  the  third  and  subsequent  terms  in  this  series  can  be 
neglected.  Substituting  equation  (12)  into  equation  (4),  we  obtain  an 
expression  of  the  form 

Z^  ■  a  +  b/y  +  additional  terms  (13) 

where  a  and  b  are  undetermined  constants,  and  the  higher  order  terms 
have  decreasing  significance. 

Because  y  is  a  function  of  Z [see  equation  (5)1,  it  is  clear 
that  the  impedance  of  Z^  placed  in  a  transmission  line  will  have  a  major 
component  that  has  the  square  root  of  the  frequency  power  dependence  of 
the  interfacial  impedance  Itself.  Thus,  for  Z^  containing  parallel 
resistance  (Rp)  and  capacitive  (Cp)  elements  only, 

1/Z.  -  1/R  +  j2itfC  (14) 

ip  p 

and  y  will  be  a  function  of  f*^;  from  equation  (13),  Zfcl  will  be  a 
function  of  f-*^,  which  is  the  form  observed. 

The  preceeding  discussion  is  Important  because  a  f” dependence 

is  often  attributed  to  diffusional  processes  at  a  plane-parallel 
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electrode.  By  mathematically  fitting  the  measured  impedance  data 

to  the  form  predicted  by  equation  (4),  we  can  clearly  see  that  the 
observed  response  is  due  to  a  capacitive  element  contained  within  a 
transmission  line,  rather  than  to  diffusion;  the  coefficients  obtained 
do  not  correspond  to  any  known  diffusional  process. 

The  capacitance  under  consideration  is  simply  that  of  the 
electrical  double  layer  on  the  sides  of  the  growing  crack.  A  resistance 
appears  in  parallel  with  this,  associated  with  Faradaic  processes  on  the 
crack  wall.  At  present,  we  are  unable  to  determine  whether  this 
Faradaic  process  is  anodic  (metal  dissolution)  or  cathodic  (reprecipita¬ 
tion  of  metal  or  hydrogen  evolution).  To  resolve  this  and  other 


important  issues,  we  intend  to  construct  a  discrete  (finite  element) 
transmission  line  model  for  the  growing  crack  based  on  the 
electromechanical  parameters  calculated  from  Impedance  measurements. 

The  two  data  sets  presented  in  Figure  8  are  not  greatly  different 
in  form;  the  major  changes  observed  are  in  the  characteristic  frequency 
and  the  intercept  with  the  real  axis  at  limiting  low  frequencies.  The 
characteristic  frequency  (fQ  -  frequency  of  minimum  phase)  moves  from 
0.14  to  0.32  Hz  between  the  first  and  second  data  sets,  and  the  limiting 
low  frequency  intercept  decreases  from  3915  to  1995  Q.  This  is  believed 
to  be  an  effect  of  crack  opening.  For  the  data  set  labeled  “crack",  the 
vestigial  crack,  although  short,  is  very  tight  and  the  impedance  is 
large.  Under  the  same  load  conditions,  the  crack  opening  displacement 
increases  as  the  crack  grows,  and  the  impedance  becomes  less.  In  fact, 
a  plot  of  limiting  resistance  versus  crack  length  shows  a  minimum  at 
about  1-cm  crack  length,  the  resistance  then  increasing  as  the  crack 
lengthens,  indicating  that  effects  of  crack  length  become  more  Important 
than  the  average  crack  opening  displacement. 

The  effects  of  crack  opening  displacement  can  be  examined  by 
varying  the  cyclic  load  conditions.  Specimens  were  cycled  at  a  mean 
load,  plus  and  minus  the  cyclic  load  perturbation.  Figure  9  shows  the 
effects  of  varying  both  the  mean  and  the  perturbation  stress  levels  for 
a  mature  crack  (*  2  cm  long)  cycled  at  0.02  Hz.  For  these  data  (and  for 
all  Impedance  data  presented  in  this  report),  each  Impedance 
determination  represents  an  average  taken  over  one  complete  period  of 
the  cyclic  load  and  thus  reflects  some  weighted  average  of  crack  opening 
displacements.  Figure  9  clearly  shows  the  effect  of  crack  opening 
(proportional  to  load)  on  the  limiting  low  frequency  resistance.  From  a 
superficial  examination,  it  would  appear  that  the  effects  of  the  load 
perturbation  level  are  more  significant  than  those  of  mean  load. 

However,  we  have  not  attempted  to  separate  the  effects  of  mean  and 
perturbation  stress  levels  quantitatively. 


Summary  and  Conclusions 


The  principal  features  of  the  AC  Impedance  response  of  a  T1-6A1-4V 
specimen  under  cyclic  load  conditions  in  3.5%  NaCl  appears  to  be  due  to 
faradaic  processes  at  the  crack  tip  and  a  transmission  line  impedance 
associated  with  the  crack  walls.  We  have  observed  qualitatively  the 
effects  of  crack  length  and  crack  opening  displacement. 


V  SECOND  YEAR  RESEARCH  PLANS 


Within  the  framework  of  the  three  research  tasks  listed  in  Section 
IV,  we  plan  to  perform  the  following  specific  subtasks  during  the  second 
year. 

(1)  Under  constant  load  conditions  (the  SCC  situation),  continue  to 
measure  AC  impedance  spectra  Inside  static  and  growing  cracks 
as  a  function  of  crack  length,  stress  intensity,  and  specimen 
potential.  Material/environment  combinations  other  than  Ti- 
6A1-4V  and  HY80  in  3. 52  NaCl  solutions  will  be  examined  to 
obtain  AC  impedance  data  for  a  range  of  SCC  behavior. 

(2)  Under  cyclic  loading  (or  CF)  conditions,  continue  to  measure  AC 
impedance  spectra  inside  static  and  growing  cracks  in  T1-6A1-4V 
and  HY80  specimens  in  3.52  NaCl  as  a  function  of  crack  length, 
AK,  R  ratio,  loading  frequency,  and  specimen  potential  under 
cyclic  loading  (or  CF)  conditions.  By  varying  AK  and  the  R 
ratio,  we  hope  to  see  the  effect  of  crack  closure  on  the  AC 
impedance  spectra.  The  electrochemical/mechanical  impedance 
will  be  measured  and  evaluated  to  see  if  it  can  help  in 
determining  the  crack  growth  mechanism.  Additional 
material/environment  combinations  will  be  examined  if  time 
permits. 

(3)  Begin  to  deconvolve  the  AC  Impedance  spectra  into 
electrochemical  parameters  that  are  pertinent  to  environment- 
assisted  cracking.  Evaluate  various  cracking  models  based  on 
this  information. 

During  Year  3  we  will  do  only  a  few  key  experiments;  most  of  our 
effort  will  be  to  use  the  data  gathered  in  Year  1  and  Year  2  to  explain 
and  model  the  electrochemical  and  mechanical  processcj  that  occur  in  SCC 
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Appendix  A 

IMPEDANCE  DATA  MANAGEMENT  AND  CONTROL  PROGRAM 


1  RER  'DNSS. RAIN  15/JUN/1983 

3  REN  BY  R.C.H.RCKUBRE 

10  HIRER:  36300:  LOREN:  24577 

90  ZNX  =  120:  IF  NPX  <  >0  THEN  1000 

100  DIN  ZR(128),ZN(5,ZNX!,ZZ(9),ZT(32).ZD*<8> 

110  PRINT  CHR4  (4) ; "BLOAD  B.FRA" 

140  ZH(1)  =  1: ZN(2>  =  1:ZR(11)  =  100:ZR112)  =  1:ZR(131  =  1:ZR(21)  =  1:ZR(22)  =  1:ZR(23)  =  1:ZR(24)  =  1: 
ZHI25)  =  1:ZN(34)  =  65535:ZR(35)  =  1:ZR(38)  =  ZNI:ZR<39>  =  1.1:ZR(44)  =  1:ZR(45)  =  0201:ZN(49>  =  2 
00 

150  ZN151)  =  1:ZR(53)  =  1:ZR(55)  =  1:ZR(56)  =  ZNX:ZR<57)  *  1:ZH(B4>  =  2:ZR(85)  =  2:ZN(86)  =  2:ZR(121)  = 
2: ZH < 122)  *  1:ZR(123)  »  1 


160  X  »  FRE  (0):SBX  *  6: 0B1  =  1:VBI  =  23 

170  IF  SI  <  4  OR  SI  >  6  THEN  SI  =  5:DX  =  1:VX  =  0 

200  60SUB  20000 

300  DATA  ■FRV'ARS’BI'/NF’.'RO’.ISV/LE'.'SR'/SL'/AC’/AS^.’AV/AE'/AP'.’OA^j'BB'/QN't'SQV’S 
S’.'RG* 

320  DATA  " IS". * IC", “RS“ , "RC", "HA*, "All", "RA" , “BE", 'DC*, “IP", “SA", "SI*, *RE“, “RA“, "RI“» "GD*, "60" , "GS" , "BR 
VNHVNSVSSVHSVSCVSOVCO' 

340  DATA  ■VE“,"ftA-,-PD","XI","XH,‘,“XL“,"XOm,"Yl","tH*,"YL","YO",-SP","PL","SH","ST","PG","VA-,-UP","D 

nvfdvfovfcvll* 

360  DATA  aBKa,atLa,aIEa,atCa,atFa,a$8a,atDa,atIa,at8a,,tPa,aEPa,aJPl,,aCPa,aTRa,aTTa,a0Pa,aFNa 
380  DATA  'AF’/BF'/RF'/TF'/UA'^UB'.'UR'/UTS'OB'.'LA'.'LB’.'LR’.’Lr.’LS'.'SV/PPVPSS^S'.’OT 
a,aRRa,aRHa 

400  DATA  "NR", ’HA* ,  PN*, ’ER", "AI" , aARa, aTSa 
1000  TEXT  :  POKE  34,0:  HOHE  :  PRINT 
1100  PRINT  :  PRINT  ‘PR06RAR  OPTIONS  ARE:" 

1120  PRINT  :  PRINT  •  1...SET  UP  SOLARTRON* 

1140  PRINT  :  PRINT  ■  2... HELP’ 

1160  PRINT  :  PRINT  a  3... TAKE  REASUREHENTS* 

1180  PRINT  :  PRINT  ■  4... CHAIN  BODE  PLOT* 

1200  PRINT  :  PRINT  "  5... CHAIN  NY8UIST  PLOT' 

1220  PRINT  :  PRINT  a  6... SAVE  DATA  ON  DISK' 

1230  PRINT  :  PRINT  a  7... RETRIEVE  DATA  FROM  DISK' 

1240  PRINT  :  PRINT  a  8... SET  TIRE* 

1245  PRINT  :  PRINT  a  9... CHAIN  ANALYSIS  PR06RAR' 

1250  PRINT  :  PRINT  a  O...EXIT  PR08RAR" 

1260  ZZ (3)  *  0:  VTAB  <23):ZZ(4)  =  0:ZZ(5)  =  9:  80SUB  10100 
1280  IF  ZEX(O)  *  1  THEN  1000 

1300  ON  ZZ (3)  GOTO  1500,2500,3000,4000,5000,6000,7000,8000.9700 

1320  HONE  :  VTAB  9:  PRINT  ‘EXIT  PROGRAM":  PRINT  :  PRINT  "TYPE  GOTO  1000  TO  RESUME":  END 

1500  60SUB  15000:  60T0  1000 

1900  ZH1128)  *  1:W  *  "T1F2":  GOSUB  29000:ZZ<4)  *  ZZ13I:RI  *  "T1F0":  60SUB  29000:  PRINT  :  PRINT  aRS  =  a 


;ZZ<4):  PRINT  :  INPUT  "USE  THIS  VALUE  (Y/N)  ?  IF  H*  <  >  *Y'  THEN  15000 
1920  ZH151)  =  ZZ(4)s  B0T0  1500s  REN  RP 
2000  HONE  :  PRINT  'INITIALISE1 
2020  60SUB  20000: H*  =  ,:,TT2'':  GOSUB  21000 
2040  60SUB  29000:  REN  KEITHLEY 
2060  GOTO  15000 

2500  HONE  :  VTAB  5:  PRINT  '  HELP:':  PRINT  *  - VTAB  10:  PRINT  '  THERE  IS  NO  HELP  AVAILABLE  HER 

E.':  PRINT  'LIFE  IS  LIKE  THAT.  IT  CAN  SOMETIMES  BE':  PRINT  'HARD  AND  CRUEL.':  PRINT  '  IF  YOU  DON 
’T  KNOH  NHAT  YOU  ARE  D0IN6' 

2510  PRINT  'IT  IS  POSSIBLE  THAT  YOU  SHOULD  NOT':  PRINT  'BE  DOING  IT.':  PRINT  •  PLEASE  DO  NOT  HESITAT 
E  TO  NOT':  PRINT  'CALL  HE,':  VTAB  19:  HTAB  9:  PRINT  'HIKE' 

2520  VTAB  23:  BET  ZD* (01:  PRINT  ZD*(0):  60T0  1000 

3000  HONE  :  VTAB  3:  PRINT  'NANE  AND  LOCATION  OF  SAVE  DATA  FILE* 

3010  60SUB  11500:  PRINT  :  PRINT  ‘FILE  NANE  NILL  BE  PRECEDED  BY  DNSS(D)-':  PRINT  :  INPUT  'ENTER  FILE  NA 
HE  =  DHSS(D)-':ZDI(0) 

3015  X  =  FRE  <0):ZF*  =  ZD*(0):ZZ<8)  =  77:ZD*(0>  =  *DHSS(P!-'  ♦  ZF*:  GOSUB  11400:ZF*  =  'DHSS(D)-'  ♦  ZF* 

+  *I':ZZ(8)  =  78 

3020  PRINT  :  PRINT  'STARTING  FILE  t  ';:ZZ( 3)  *  ZFX:ZZ<4)  =  0:ZZ(5)  =  100:ZZ(6)  =  1:  60SUB  10100:  IF  ZE 
X(0)  =  1  THEN  3020 
3030  ZFX  =  ZZ(3):ZN(3)  =  ZHI1181 

3040  X  *  FRE  (0):  POKE  34,0:  HONE  :  POKE  34,3:  VTAB  1:  PRINT  'NEASUREHENT  NODE  »  ’;ZN(127):  VTAB  2:  PRINT 
ZF*; ZFX:  VTAB  3:  PRINT  'SCALE  FACTOR  ':ZN(53) 

3050  ZHI47)  =  0:  HONE  :  VTAB  5:  PRINT  'STATUS  PAGE':  PRINT  * . ':  VTAB  8:  PRINT  'TOTAL  I  OF  POI 

NTS  ':ZNX:  VTAB  10:  PRINT  'TIHED  SWEEP';:  HTAB  19:  PRINT  'ON':  IF  ZH1125I  =  0  THEN  VTAB  10:  HTAB 
19:  PRINT  'OFF* 

3060  VTAB  12:  PRINT  'DC  VOLTAGE  SWEEP';:  HTAB  19:  PRINT  'ON*:  IF  ZHI120)  =  0  THEN  VTAB  12:  HTAB  19:  PRINT 
•OFF' 

3070  VTAB  14:  PRINT  'FREQUENCY  SHEEP'; :  HTAB  19:  PRINT  'ON':  IF  ZN134)  *  ZHC35)  THEN  VTAB  14:  HTAB  19 
:  PRINT  'OFF' 

3080  VTAB  16:  PRINT  'HARHONIC  SHEEP';:  HTAB  19:  PRINT  'ON':  IF  ZNI122)  =  ZHU23)  THEN  VTAB  16:  HTAB  1 
9:  PRINT  'OFF* 

3090  VTAB  18:  PRINT  "INPUT  FRONT/REAR';:  HTAB  19:  PRINT  'FRONT':  IF  ZHI50)  >  0  THEN  VTAB  18:  HTAB  19: 

PRINT  'REAR  ':  IF  ZH150)  >  1  THEN  VTAB  18:  HTAB  19:  PRINT  'BOTH' 

3100  VTAB  2:  HTAB  28:  PRINT  'START  9  -;ZM(124):  VTAB  3:  HTAB  29:  PRINT  'STOP  9  ';ZH(125> 

3110  IF  ZHI126)  >  0  THEN  GOSUB  18200:  IF  ZZ(3>  ♦  ZZ(4>  /  60  <  ZNI124)  THEN  FOR  II  =  1  TO  4444:  NEXT 
II:  60SUB  18000:  60T0  3110 
3120  FRI  =  1:  IF  ZHI50)  =  2  THEN  FRX  =  2 
3130  ZN (56)  =  INT  (ZNX  /  (FRX  I  ZH<55)  I  ZH<57))) 

3140  ZH (39)  =  (ZH (34)  /  ZH(35H  A  (1  /  ZH(56M 

3145  ZN(O.O)  =  ZNX:  GOSUB  18200: ZN(l.O)  *  ZZ(3)  ♦  ZZ(4!  /  60:ZN(2,0)  =  ZM4127) : ZN(3,0)  =  ZH(53):ZN(4,0) 

=  1250 

3150  VTAB  23:  PRINT  ZH(39),ZH(57),ZH<56> 

3160  VTAB  5:  CALL  -  958:  PRINT 

3170  J  -  0:  IF  ZH ( 121 )  *  1  THEN  ZH(3>  =  ZH(llB):  REN  INITIAL  DC 

3175  GOSUB  3900:  REN  SET  DISPLAY 

3180  II  =  3:  GOSUB  17000 

3190  FOR  JZN  *  1  TO  441  I  ZHU17):  NEXT  JZN 

3200  ZH(1)  *  ZH (35) :  REN  F  NIN. 

3210  II  -  1:  60SUB  17000 

3220  ZH129)  *  0:  IF  ZH(1>  >  =  ZHI49)  THEN  ZN(29)  =  1:  REN  COUPLING 

3230  II  *  29:  GOSUB  17100 

3240  ZH(2S)  =  ZH(122):  REN  INITIAL  HARHONIC 


3230  II  =  25:  60SUB  17000 

3260  FRZ  =  0: ZH (30)  =  Os  IF  ZHI50)  =  1  THEN  ZH(30)  =  IsFRZ  =  Is  REN  REAR 

3270  II  =  30s  60SUB  17100:ZH(45)  =  0100s  IF  ZHI127)  =  2  THEN  ZH(45)  =  0200 

3280  IF  ZN1127)  =  0  THEN  ZH(4S)  =  0201 

3280  J  *  3  ♦  Is  IF  J  =  ZNZ  AND  ZH<128)  *  1  THEN  3620 

3283  IF  J  >  ZNZ  THEN  3380s  REN  FULL 

3300  II  =  45:  6GSUB  17000s  60SUB  18000s II  =  32:  6DSUB  16800s  REN  HEASURE 
3310  HI  =  *:,0P2,r*sII  =  4s  60SUB  24000s  IF  ZZ<4)  <  >  0  THEN  3300 

3320  FOR  ZS  =  1  TO  3:ZN(ZS  -  1,J)  =  ZZ(ZS):  NEXT  ZS 

3330  IF  ZNI127)  >  3  THEN  II  =  4S:ZH(45)  =  0200:  60SUB  17000:  60SUB  IBOOOsNI  =  V0P2,r'iII  =  4:  60SUB 
24000:  REH  CHI2 

3340  IF  ZN(127)  <  3  THEN  ZN(3,J)  =  1:ZN(4,J)  =  0 

3330  IF  ZNI127)  =  3  THEN  ZN(3,J)  *  ZN ( 1 , J) : ZN (4, J >  =  ZN(2,J):ZN(1,J)  =  ZZ(2):ZN(2,J)  =  ZZ(3) 

3360  IF  ZN ( 127)  *  4  THEN  ZN(3,J)  *  ZZ(2)  -  ZN(1,JI:ZN(4,J)  =  ZZ( 3)  -  ZN(2,J) 

3370  IF  ZHI127)  =  3  THEN  ZN(3,J)  =  ZN( 1 , J ) : ZN ( 4 , J )  -  ZN (2, J> s ZN<1 , J)  =  ZZ (2) : ZN (2, J)  *  ZZ(3) 

3380  IF  ZNI127)  =  5  THEN  60SUB  3800:  60T0  3420 

3380  60SUB  9500:ZN(3,J)  =  1  /  ZN(51) : ZN<4, J)  =  6.28232  t  ZN(0,J)  I  ZN452) :  80SUB  9500:ZN(1,J)  *  ZN(1,J 
)  /  ZN(53):ZN(2,J)  *  ZN(2,J)  /  ZH(S3)s  REH  SCALE 
3400  GOSUB  9980:ZN(3, J)  *  ZH(3):  REH  HAG . IrPHASE 
3420  II  =  1:H$  =  ■s*?ER,,s  60SUB  24000:  VTAB  21:  HTAB  12:  PRINT  ZZ(II) 

3430  ZZ (3)  =  7: ZZ <4>  =  ZH(  122) : ZZ (5)  =  ZH(U8) :ZZ<6)  *  ZH(35):  GOSUB  3990 
3432  ZZ(3)  =  8: ZZ (4)  =  ZH( 123) : ZZ (5)  *  ZH< 1 19) : ZZ (6)  *  ZH134):  GOSUB  3990 
3434  ZZ(3)  *  9s ZZ (4)  =  1:ZZ(5)  =  ZH ( 1 20) : ZZ (6)  =  ZH139):  60SUB  3990 
3436  ZZ (3)  *  10: ZZ (4)  =  ZH(25):ZZ(5)  *  ZH(3):ZZ(6)  =  ZH(i):  GOSUB  3990 
3440  VTAB  12:  HTAB  9:  PRINT  J  -  IjVjZNX;  TAB<  25) ; J; ■/•; ZNZ 

3450  FOR  II  =  0  TO  5:  VTAB  (13  ♦  II):  HTAB  6:  CALL  -  868:  PRINT  ZNCII.J  -  1);  TAB (  24) ; ZN( II . J) :  NEXT 
II 

3460  INVERSE  :  VTAB  20:  HTAB  !2:ZDI<0)  =  'FRONT*:  IF  ZH(30)  *  1  THEN  ZDKO)  =  *  REAR' 

3462  PRINT  ZDKO);:  HTAB  28:ZDK0>  *  'DC':  IF  ZH(29)  *  1  THEN  ZDKO)  *  'AC' 

3464  PRINT  ZDKO);:  HTAB  36:ZDK0)  =  'OFF':  IF  ZH(14)  =  1  THEN  ZDKO)  =  'ON  * 

3466  PRINT  ZDKO):  NORHAL 

3480  VTAB  4:  CALL  -  868:  IF  ZH(47)  *  999  THEN  PRINT  'TERHINATE' 

3490  IF  ZH(50)  =  2  AND  FR1  *  0  THEN  ZHI30)  =  1:FRI  =  1:  60T0  3270:  REH  REAR 
3500  IF  ZH(25)  <  ZH(123)  THEN  ZH(25)  =  ZH(25)  ♦  1:  GOTO  3250:  REH  LOOP  H 
3510  IF  ZH<1)  >  =  ZH(34)  THEN  3540 

3520  ZHI1)  =  ZH(l)  I  ZH(39):  IF  ZH(1)  >  ZH(34)  THEN  ZH(1)  =  ZHC34) 

3530  60T0  3210 

3540  IF  SGN  <ZH(120>)  «  S6N  (ZHU19)  -  ZH(3))  THEN  ZH(3)  *  ZH(3)  *  ZHII20):  GOTO  3180 

3530  REH  DC  LIHIT  REACHED 

3560  IF  ZH(121)  -  1  THEN  3180 

3570  GOTO  3210 

3580  REH  FULL 

3590  GOSUB  18000 

3600  ZH ( 124)  *  ZHU24)  *  ZH<126>:  REH  INC.TIHE 

3610  X  *  FRE  (0):ZFX  *  ZFX  ♦  lsZDKO)  *  ZFt  ♦  STRI  (ZFZ):ZZ(8)  *  78:  GOSUB  11400:  REH  SAVE  DATA 
3612  IF  <ZZ (3)  ♦  ZZ (4)  /  100)  >  ZHI125)  OR  Zfl<47)  -  999  THEN  1000:  REH  TERHINATE 
3616  IF  ZH<121)  *  2  THEN  ZH(3I  -  ZH<3)  ♦  ZH(120):  IF  S6N  (ZH(12))  •  S6N  (ZH(U9)  -  ZH(3) )  THEN  1000: 
REH  DC  LIHIT 

3618  80T0  3040:  REH  CLOSE  LOOP 

3620  ZH(1)  «  60:11  ‘  1:  GOSUB  17000: ZNIO, ZNZ)  *  60s  FOR  A  *  1  TO  441:  NEXT  A:Hf  *  "TIFl's  GOSUB  29200:Z 
N12.ZNZ)  -  ZZ(3):HI  *  'TlFO's  GOSUB  29200: ZN(3, ZNZ)  «  ZZ(3) 


3630  ZN ( 18)  =  0:!1  =  18:  80SUB  17000:  FOR  II  =  1  TO  441:  NEXT  II:H*  =  "TlFl “ :  80SUB  29200: ZN(4, ZNZ)  =  2 
Z(3):NI  =  •T1F0":  60SUB  29200: ZN(5,ZNX)  =  ZZI31 
3640  II  =  20:  60SUB  16900:  80T0  3420 
3800  REN  CH1=I,CH2=V 

3820  ZNI3.J)  =  ZN (25) : ZN (4, J )  =  ZZ(2):ZN(5,J)  =  ZZI3) 

3840  ZZ (2)  •  1  /  ZN(51> :ZZ<3)  =  6.28232  t  ZN(O.J)  t  ZH(52):ZZ(4)  =  ZNIl.Jl  t  ZZI2)  -  ZN(2,J)  t  ZZ(3I:ZN 
(2, J)  =  ZN(l.J)  t  ZZ (3)  ♦  ZN (2, J )  I  ZZ (2) : ZN < 1 , J )  =  ZZI4):  REN  I=E/R 
3860  ZN(1,J)  =  ZN ( 1 , J )  t  ZN(53) : ZN(2, J)  =  ZN(2,J)  t  ZN(53) : ZN (4, J)  =  ZN(4,J)  t  ZN (53) : ZN (5. J )  =  ZNI5.J) 
I  ZN(53> :  REN  SCALE  NULTIPLIER 
3880  RETURN 
3900  REN  DISPLAY 
3920  REN 

3930  VTAB  6:  CALL  -  958:  HTAB  8:  PRINT  *H  DCV  FREQUENCY*:  PRINT  *INIT.*:  PRINT  'FINAL*:  PRINT  * 
STEP*:  PRINT  *CURR.*:  FOR  II  =  0  TO  5:  VTAB  (13  ♦  II):  PRINT  *#*;II:  NEXT  II 
3940  VTAB  20:  PRINT  ‘INPUT  FRON  :  C0UPLIN6  :CONP  * 

3960  VTAB  21:  PRINT  'LAST  ERROR:  *:  RETURN 

3990  VTAB  ZZ(3):  HTAB  7:  CALL  -  868:  PRINT  ZZ<4);  TAB <  tl);ZZ(5);  TAB (  18);ZZ(6>:  RETURN 
4000  ZDKO)  =  'BODE*:  60SUB  19900 
4020  NPX  =  1:  CALL  520*DNSS.B0DE* 

5000  ZD$(0)  *  'NYQUIST':  60SU8  19900 
5020  NPX  =  1:  CALL  520'DNSS.NYQUIST* 

6000  REN  SAVE  DATA  TO  DISK 

6010  60SUB  11500:  PRINT  :  PRINT  ‘FILE  NANE  NILL  BE  PRECEDED  BY  DHSS(D)-*:  PRINT  :  INPUT  'ENTER  FILE  NA 
HE  =  DHSS(D)'*;ZD$(0):ZDt(0)  =  'DNSS(D)-*  *  ZD*(0) 

6020  ZZ (8)  =  78:  60SUB  11420:  6OT0  1000 

6500  ZDt(O)  =  ZD* (0)  ♦  STR*  (  INT  (ZZ(3)  ♦  0.5)):  RETURN 

7000  REN  RETRIEVE  DATA  FRON  DISK 

7040  X  =  FRE  (0):  60SUB  11500:  PRINT  :  PRINT  ’CURRENT  FILE  NANE  «  *;ZF»;*»*;ZF1:  PRINT  :  INPUT  *  NEH 
FILE  NANE  =  DHSS(D)-*;ZD*(0):ZDII0>  =  *DNSS(DI-*  ♦  ZDKO):  IF  ZDt(O)  <  >  *DNSS(D)-'  THEN  ZF*  =  Z 
D*  (0) 

7060  X  =  FRF  (0):  INPUT  *  NEH  FILE  •  =  *;ZFZ:ZD«(0)  *  ZF*  ♦  *»*  ♦  STR*  (ZFI) 

7080  ZZ(B)  *  78:  60SUB  11300:  60T0  1000 

8000  HONE  :  REN  SET  TINE  _ 

8020  60SUB  18200:  VTAB  8:  PRINT  'NEH  TINE:*:  HTAB  9:  INPUT  'HOURS  *;JZN:  IF  JZN  >  99  OR  JZN  <  0  THEN  8 

020 

8040  VTAB  11:  HTAB  7:  INPUT  'NINUTES  *:ZN(83):  IF  ZR(83)  )  59  OR  ZH(B3)  <  0  THEN  8040 
8060  II  -  83:  60SUB  17200:  60T0  1000 
9000  REN  DISPLAY  COUNTER  ROUTINE 

Qfilfl  UTAH  TA*  PAINT  *t  '•J7N>a  *•  RFTURN 

9240  PRI  3:  PRINT  *8(s*t  PRINT  *’*;ZD»(0);*X'*:  PRINT  *i?;8H:*:  PRI  0:  INPUT  ZD4(0):  IM  0:  RETURN 

9500  REN  COHPLEXDIVISION 

9520  ZZ(5)  <  ZN (3, J )  »  ZN(3,J)  4  ZN(4,J)  t  ZNI4.J) 

9540  ZZ (6)  «  ( ZN ( 1 « J )  I  ZN(3,J)  4  ZN(2,J)  t  ZN(4,J)>  /  ZZ(5) 

9560  ZZ (4)  *  (ZN(2,J)  I  ZN(3,J)  -  ZN(1,J)  t  ZN(4,J))  /  ZZ(5) 

9580  ZN(1,J)  >  ZZ(6):ZN(2,J)  =  ZZ(4):  RETURN 
9700  ZDKO)  *  'ANALYSIS*:  60SUB  19900 
9720  CALL  520'DNSS.ANALYSIS* 

9970  REN  NAS.  6  PHASE 

9980  ZN(3,J)  >  SQR  (ZNd.J)  I  ZN(1,J)  4  ZN(2,J)  t  ZN(2,J)) 

9985  ZN(4,J)  *  ATN  (ZN(2,J)  I  ZN(l.J))  I  <180  /  3.141592654) 

9990  ZN(5,J)  *  L06  (ZNI3.JH  /  LOS  (10):  RETURN 
9999  PRINT  CHR*  (13);  CHRI  (4);*SAVEDNSS.HA1N*:  END 


10000  REN  Y/N  SUBR 

10010  ZENO)  =  0:  PRINT  *VES  ";:  BET  ZDI(0):ZZ(3>  =  PEEK  (36)  -  4:  IF  ZDKO)  =  "Y"  OR  ZDKO)  =  CHRt  ( 

13)  THEN  ZZ(3)  =  2:  PRINT  "":  RETURN 

10020  IF  ZDKO)  =  "N"  THEN  POKE  36.ZZ13):  PRINT  ’NO  -sZZ13)  =  Is  RETURN 

10030  ZENO)  =  1:  60SUB  10050:  RETURN 

10050  REN  ERROR  SUBR 

10060  PRINT  CHRI  (7):  ZENO)  =  1:  RETURN 

10100  REN  I/P  SUBR 

10110  CALL  -  868 

10120  ZENO)  =  0:  PRINT  ZZ< 3):  INPUT  *  a;Z0*(0):  IF  ZDKO)  =  ""  THEN  ZZ(7)  =  ZZ( 3):  BOTO  10170 
10130  FOR  I  -  1  TO  LEN  (ZDKOOsZBI  *  NIDI  (ZDKO), 1,1) 

10140  IF  ZOI  =  OR  ZOI  *  *-■  OR  ZBI  =  *  *  OR  ZOI  =  V  OR  ZD!  =  "E"  THEN  10160 

10150  IF  ZOI  <  *0*  OR  ZOI  >  "9"  THEN  10050 

10160  NEXT  I: ZZ  (7)  =  VAL  (ZDKO)) 

10170  IF  ZZ (7)  <  ZZ (4)  OR  ZZ< 7)  >  ZZ<5)  THEN  10050 
10180  ZZ(3)  *  ZZ (7) :  IF  ZZ(6)  *  1  THEN  ZZ(3)  *  INT  CZZ<7>) 

10190  FOR  I  =  0  TO  200:  NEXT  I:  RETURN 

11300  REN  :  LOAD  ARRAY 

11320  POKE  98,90:  POKE  99,ZZ(8):  CALL  38152 

11330  PRINT  CHRt  (4);'BL0ADaZDI(0);',Aa:  PEEK  (38394)  ♦  PEEK  (38395)  I  256 

11340  RETURN 

11400  REN  :  SAVE  ARRAY 

11420  POKE  98,90:  POKE  99,ZZ(8):  CALL  38104 

11430  PRINT  CHRI  (4) ; "BSAVt'ZDI (0) ",  A"  PEEK  (38394)  ♦  PEEK  (38395)  I  256', L*  PEEK  (38396)  ♦  PEEK  (3 
8397)  I  256* 

11440  RETURN 

11500  REN  ASK  SLOT, DRIVE  AND  VOLUHE 

11510  VTAB  6:  CALL  -  958:  PRINT  ’POINTERS  CURRENTLY  SET  AT  :a:  VTAB  8:  PRINT  a  SLOT  #a;SZ:  PRINT  a 
DRIVE  #a;DX:  PRINT  a  VOLUHE  f";VZ 

11520  VTAB  12:  PRINT  ‘CHAN6E  POINTERS  (Y/N)  ?  "j:  8ET  ZDNO):  IF  ZDNO)  =  aNa  THEN  11550 
11530  PRINT  :  PRINT  :ZZ(3)  *  SZ:ZZ(4)  =  4:ZZ(5)  -  6:ZZ(6)  *  0:  PRINT  '  SLOT  60SUB  10120:  IF  ZEX 

(0)  -  1  THEN  11530 
11532  SX  =  INT  (ZZ (3)  ♦  0.5) 

11535  PRINT  : ZZ (3)  =  DX:ZZ(4>  =  1:ZZ(5>  =  2:ZZ(6)  =  0:  PRINT  a  DRIVE  6DSUB  10120:  IF  ZEX(O)  =  1  THEN 
11535 

11537  DZ  =  INT  (ZZ(3)  ♦  0.5):  IF  SX  <  >6  THEN  VX  =  0:  BOTO  11550:  REN  FLOPPY 
11540  PRINT  :ZZ(3)  *  V1:ZZ(4>  *  1:ZZ(5)  =  67:ZZ(6)  *  0:  PRINT  a  VOLUHE  60SUB  10120:  IF  ZEX(O)  *  1 
THEN  11540 

11542  VZ  -  INT  (ZZ (3)  *  0.5) 

11550  PRINT  :  PRINT  "PRESS  ANY  KEY  TO  CATAL06  a:  6ET  ZDKO):  PRINT  :  PRINT  CHRI  <4>;aCATAL06,SaSZa,Da 
DXa,VaVI 

11560  PRINT  :  PRINT  "ANOTHER  CATAL06  Y/N  ?  60SUB  10000:  IF  ZEX(O)  =  1  THEN  11560 
11570  IF  ZZ (3)  *  2  THEN  11500 
11580  RETURN 

15000  TEXT  :  POKE  34,0:  HONE  :  VTAB  5:  PRINT  "SET-OP  OPTIONS  ARE:" 

15020  PRINT  :  PRINT  "  1... RETURN  TO  RAIN  HENU" 

15040  PRINT  :  PRINT  "  2...CREATE/H0DIFY  SET-UP  FILE" 

15060  PRINT  :  PRINT  "  3... RETRIEVE  SET-UP  FILE  FROH  DISK" 

15080  PRINT  :  PRINT  "  4... SAVE  SET-UP  FILE  TO  BISK" 

15100  PRINT  :  PRINT  ■  5... EXECUTE  SET-UP  FILE" 

15120  PRINT  :  PRINT  "  6... RESET  (INITIALIZE)  SOLARTRON" 

15140  PRINT  :  PRINT  "  7.. .RESET  (INITIALIZE)  KEITHLEY* 


15200  ZZI3)  =  1:  VTAB  <23):ZZ<4)  =  1:ZZ(S)  =  7:  60SUB  10100:  IF  ZEZ(O)  =  1  THEM  15000 
15220  ON  ZZ (3)  GOTO  1000, 15300, 16600, 16640, 16700,2000, 1900 
15300  HONE  :  VTAB  5:  PRINT  'SET-UP  WHICH  SUBSYSTEM: ' 

15320  PRINT  :  PRINT  '  1... GENERATOR* 

15340  PRINT  :  PRINT  '  2... ANALYSER* 

15300  PRINT  :  PRINT  '  3. . .SWEEP(S) * 

15380  PRINT  :  PRINT  '  4. . . IMPEDANCE  STANDARD* 

15400  PRINT  :  PRINT  '  5... MEASUREMENT  CONFIGURATION* 

15440  PRINT  :  PRINT  *  0... RETURN  TO  SET-UP  OPTIONS* 

15480  ZZ< 3)  *  0:  VTAB  <23!:ZZ(4)  =  0:ZZ(5)  =  5:  60SUB  10100:  IF  ZEZ(O)  =  1  THEN  15300 
15500  ON  ZZ (3)  60T0  15400,15800,14000,14300,14500 
15520  60T0  15000 

15400  HOME  :  VTAB  3:  PRINT  “GENERATOR" 

15410  VTAB  4:  PRINT  -FREQUENCY  ';:ZZ<3>  =  ZM(1):ZZ<4>  =  IE  -  5:ZZ<5)  =  45535:ZZ(4)  =  0:  60SUB  10100:  IF 
ZEX(O)  *  1  THEN  15410 
15420  ZH(1)  =  ZZ (3> 

15430  VTAB  9:  PRINT  'AMPLITUDE  *;:ZZ<3)  *  ZM(2):ZZ(4)  =  0: ZZ(5)  -  10.23:ZZ(4)  =  0:  60SUB  10100:  IF  ZEX 
(0)  =  1  THEN  15430 
15440  ZH(2)  -  ZZ (3) 

15450  VTAB  12:  PRINT  *D.C.  BIAS  *;:ZZ(3)  -  ZM(3):ZZ(4)  -  -  10.23:ZZ(5)  *  10.23:ZZ(4)  *  0:  GOSUB  10100 
:  IF  ZEZIO)  =  1  THEN  15450 
15440  ZN(3)  -  ZZ (3) 

15444  IF  ZN(2)  ♦  ABS  (ZH(3>>  >  14  THEN  VTAB  8:  PRINT  ‘COMBINED  VOLTAGE  TOO  LAR6E*:  60T0  15430 
15470  VTAB  15:  PRINT  'NAVEFORM  (SIN=O.SflRM)  ';:ZZI3>  =  ZM(4):ZZ(4>  *  0.-ZZI5)  *  1:ZZ(4>  *  1:  GOSUB  101 
00:  IF  ZEZ(O)  *  1  THEN  15470 
15480  ZH(5>  *  ZZ (3) 

15490  VTAB  18:  PRINT  'AMPLITUDE  COMP.  (0FF=0,0N=1)  *;:ZZ(3>  -  ZM(14):ZZ(4I  »  0:ZZ(5)  =  1:ZZ(4>  *  1:  60SUB 
10100:  IF  ZEZ(O)  =  1  THEN  15490 
15700  ZH(14>  =  ZZ (3) 

15710  IF  ZH114)  =  0  THEN  15300 

15720  HOME  :  VTAB  5:  PRINT  ‘AMPLITUDE  DHPRESSION  ON* 

15730  VTAB  8:  PRINT  'SOURCE  CHANNEL  *;:ZZ(3)  *  ZM(11):ZZ<4)  -  1:ZZ(5)  -  2:ZZ(4>  *  1:  60SUB  10100:  IF  Z 
EZIO)  =  1  THEN  15730 
15740  ZNI11)  *  INT  1100  I  ZZI3)  ♦  0.5) 

15750  VTAB  12:  PRINT  'VOLTAGE  LEVEL  *;:ZZ<3)  =  ZH(12):ZZ(4»  =  IE  -  4:ZZ(5)  =  300:ZZ(4>  =  0:  GOSUB  1010 
0:  IF  ZEZIO)  =  1  THEN  15750 
15740  ZMI12)  =  ZZ13) 

15770  VTAB  14:  PRINT  'ERROR  PERCENT  V.ZZI3)  =  ZM(13):ZZI4)  =  1:ZZ(5>  =  50:ZZI4)  *  0:  GOSUB  10100:  IF 
ZEZIO)  =  1  THEN  15770 
15780  ZMI13)  =  ZZI3) 

15790  GOTO  15300 

15800  HOME  :  VTAB  3:  PRINT  'ANALYSER' 

15810  VTAB  4:  PRINT  'INTEGRATION  TINE  (SECONDS)  'j:ZZ<3>  -  ZMI21):ZZ(4>  *  ,01:ZZ15)  -  1E4:ZZI4)  *  0:  GOSUB 
10100:  IF  ZEZIO)  *  1  THEN  15810 
15820  ZMI21)  -  ZZ (3) 

15830  VTAB  9:  PRINT  'MEASUREMENT  DELAY  ':  GOSUB  15995:ZZ<3)  »  ZH(II):ZZ(4)  =  0:ZI(5)  *  1E5:ZZ(4)  *  0;  GOSUB 
10100:  IF  ZEZIO)  -  1  THEN  15830 
15840  ZMIII)  »  ZZI3) 

15850  VTAB  14:  PRINT  'INITIAL  HARMONIC  *;:ZZ(3)  =  ZM(25):ZZ<4)  *  1:ZZ<5)  =  14:ZZ(4)  «  1:  GOSUB  10100:  IF 
ZEZIO)  *  1  THEN  15850 
15840  ZMI2S)  *  ZZ (3) 


15870  VTAB  18s  PRINT  'AUTO  INTEGRATE  0FF*0‘:  PRINT  ‘CH1(LM,CH2(L>=2,CH1<S>=3,CH2(S)M‘:ZZ<3>  =  ZH<2 
6) s ZZ (4)  *  0: ZZ (5)  -  4:ZZ(6>  *  It  GOSUB  10100s  IF  ZEZIOl  *  1  THEN  15870 
15880  ZH<26>  •  ZZ (3) 

15900  VTAB  6s  CALL  -  958s  VTAB  8s  PRINT  ‘INPUT  V0LTA6E  RANGE  AUT0=0‘s  PRINT  *30mM,300HV=2,3V=3,30V 
=4,300V=5‘:ZZ(3)  *  ZN(27)sZZ(4)  =  OsZZ(S)  *  5sZZ(6)  *  It  60SUB  10100s  IF  ZEZ(O)  *  1  THEN  15900 
15910  ZHI27)  =  ZZ (3) 

15920  VTAB  13s  PRINT  ‘INPUT  C0UPLIN6  ‘s  PRINT  *DC«0,  ACM.  FREQUENCY  TO  CHAN6E«m*:ZZ(3)  =  ZXZ(l)sZZ(4) 

=  Os  ZZ  (5)  =  999:  ZZ  (6)  =  0:  80SU1  10100s  IF  ZEZ(O)  *  1  THEN  15920 
15930  ZNI49)  =  ZZ(3)sZH(29)  =  ZH(491s  IF  ZH1491  >  I  THEN  ZN1291  *  0 

15950  VTAB  18s  PRINT  'INPUT  (FR0NT«0, REARM, B0TH»2)  ‘; sZZ(3)  *  ZH(50tsZZ(4)  *  0:ZZ(5)  »  2:ZZ(6>  *  1:  60SUB 
10100s  IF  ZEZ(O)  *  1  THEN  15950 

15960  ZH150)  *  ZZ (3)  t Z71  (30)  *  ZN(50ls  IF  ZN150I  >  1  THEN  ZN(30)  =  0 
15990  60T0  15300 

15995  PRINT  ‘(S=SECONDS,C*CYCLES)  ‘ssll  =  24s  PRINT  ‘CYCLES  ‘ss  BET  ZDt(O):  IF  ZD410)  *  *S*  THEN  HTAB 
22:  VTAB  10:  PRINT  'SECONBS'ssIl  *  23 

15996  RETURN 

16000  HONE  :  VTAB  3:  PRINT  ‘SNEEP(S)...*  OF  NEA8URENENTS  »  ‘sZNZ 

16010  VTAB  6:  CALL  -  958:  PRINT  *0C  VOLTAGE  SHEEP  <Y/N>  ?  GOSUB  10000: ZH1571  *  1:  IF  ZEZ(O)  *  1  THEN 
16010 

16020  IF  ZZ (3)  »  1  THEN  ZH(118>  «  ZN(3>:ZH(119)  *  ZH(3):ZH(120)  >  0:  60TO  16100 
16030  VTAB  9:  PRINT  ‘INITIAL  DC  VOLTAGE  ‘;:ZZ(3>  *  ZH(3):ZZ<4)  «  -  10.23:ZZ<5)  *  10.23:ZZ(6)  =  0:  GOSUB 
10100:  IF  ZEZ(O)  =  1  THEN  16030 
16040  ZH4118I  *  ZZ (3> 

16050  VTAB  12:  PRINT  ‘FINAL  DC  VOLTAGE  ';:ZZ(31  *  ZN< 1 19) i ZZ (4)  =  -  10.23:ZZ(5)  *  10.23:ZZ(6)  =  0:  GOSUB 
10100:  IF  ZEZ(O)  *  1  THEN  16050 
16060  ZN1119)  *  ZZ (3) 

16070  VTAB  15:  PRINT  ‘DC  VOLTAGE  STEP  ‘;:ZZ(3>  *  Zlf ( 120) s ZZ 14)  *  -  10.23tZZ(5)  «  10.23:ZZ(6>  =  0:  60SUB 
10100:  IF  ZEZ(O)  *  1  THEN  16070 
16080  ZN1120)  =  ZZ (3) 

16085  IF  SGN  (ZN(U9)  -  ZfllllSI)  <  >  S6N  (ZH11201)  THEN  VTAB  8:  PRINT  'ARE  YOU  60IN6  TO  CIRCLE  INF 
INITY  ?*:  GOTO  16030 

16088  IF  ZH1120)  *  0  THEN  ZHU21)  *  2sZH(57)  =  lsZN4117)  -  0:  GOTO  16100 

16090  VTAB  18:  PRINT  'FILL  ARRAY  BETNEEN  DC  STEPS  (Y/N)‘;:  GOSUB  10000:ZH(121)  -  ZZ13):  IF  ZH(121)  =  2 
THEN  16100 

16093  ZH157)  *  INT  (1.5  ♦  (ZNU19)  -  ZH(U8>>  /  ZH(120>> 

16095  VTAB  20:  PRINT  ‘DELAY  BETNEEN  DCV  </S)  ‘;:ZZ(3)  *  ZH(U7):ZZ(4>  *  0:ZZ(5I  =  1000:ZZ(6)  =  0:  60SUB 
10100:  IF  ZEZ(O)  *  1  THEN  16095 
16097  ZH(117)  *  ZZ(3> 

16100  VTAB  6:  CALL  -  958:  PRINT  'FREQUENCY  SHEEP  (Y/N)  ?  60SUB  10000:ZN(56)  =  1:  IF  ZEZ(O)  «  1  THEN 
16100 

16105  IF  ZZ (3)  =  1  THEN  ZH(341  =  ZH(l)slI  -  34:  60SUB  21000:ZH(35)  =  ZNIllsII  =  35:  GOSUB  21000:  60T0 
16150 

16110  VTAB  9:  PRINT  ‘FREQUENCY  HINUHUH  *;:ZZ(3)  •  ZN(35):ZZ(4)  =  IE  -  5:ZZ<5)  *  65535:ZZ(6)  =  0:  60SUB 
10100:  IF  ZEZ(O)  -  1  THEN  16110 
16120  ZH43S)  »  ZZ(3):II  *  35:  GOSUB  17000 

16130  VTAB  12:  PRINT  'FREQUENCY  NAZINUH  ';:ZZ(3)  -  Z99434) : ZZ (4)  -  ZH(35):ZZ(5)  =  65535:ZZ(6)  =  0:  GOSUB 
10100:  IF  ZEZ(O)  -  1  THEN  16130 
16140  ZH<34)  -  ZZ(3):II  •  34:  60SUB  17000 

16150  ZN(55)  =  1:  IF  ZN(34)  >  32767  THEN  ZN(25)  -  1:  GOTO  16200:  REN  NO  HARHONIC  POSSIBLE 
16155  VTAB  6:  CALL  -  958:  PRINT  ‘HARHONIC  SHEEP  (Y/N)  ?  GOSUB  10000:  IF  ZEZ(O)  =  1  THEN  16150 
16160  IF  ZZ(3)  «  1  THEN  ZHU22)  =  ZH(25):ZN<123>  »  ZH(25):  GOTO  16200 


16170  VTAB  9:  PRINT  'INITIAL  HARMONIC  I  ";:ZZ(3>  =  ZH ( 122) : ZZ (4)  =  1:ZZ(S)  =  INT  (6SS3S  /  ZH(34)):ZZ( 

6)  =  1:  GOSUB  10100:  IF  ZEKO)  =  1  THEN  16170 
16180  ZH1122)  *  ZZ (3) 

16190  VTAB  12:  PRINT  'FINAL  HARMONIC  •  *;:ZZ( 3)  =  ZM< 123) : ZZ (41  =  ZM ( 122) : ZZ (5)  =  INT  (65535  /  ZM134) 

) :  ZZ  (6)  =  1:  60SUB  10100:  IF  ZEKO)  =  1  THEN  16190 
16195  ZM( 123)  =  ZZ (3) : ZM(55)  =  ZH<123>  -  ZM(122)  ♦  1 

16200  VTAB  6:  CALL  -  958:  PRINT  -TIMED  SNEEP  (Y/N)  ?  ';:  60SUB  10000:  IF  ZEKO)  =  1  THEN  16200 
16210  IF  ZZ (3)  =  1  THEN  ZHI124)  =  0:ZM(125)  =  ZM(124) : ZMI126)  =  ZN(125):  GOTO  16290 
16220  60SUB  18000:  VTAB  9:  PRINT  ZDKO):  REN  CLOCK 

16230  VTAB  12:  PRINT  ‘START  MEASUREMENTS  9  HH.HH  ‘;:ZZ(3>  =  ZMU24) : ZZ (4)  *  0:11(5)  =  24:ZZ(6)  *  0:  GOSUB 
10100:  IF  ZEKO)  =  1  THEN  16230 
16240  ZH( 124)  =  11(5) 

16250  VTAB  15:  PRINT  'STOP  MEASUREMENTS  9  HH.HH  *;:ZZ(3)  =  ZH( 125) : ZZ (4)  =  ZH( 124) : ZZ (5)  *  99.59:ZZ(6I 
=  0:  GOSUB  10100:  IF  ZEKO)  =  1  THEN  16250 
16260  ZH(125)  =  11(5) 

16270  VTAB  18:  PRINT  ‘NEASUREHENT  INTERVAL  HH.HH  '*:ZZ< 3)  =  ZH ( 126) : ZZ (4)  =  0:11(5)  =  24:ZZ(6)  *  0:  60SUB 
10100:  IF  ZEKO)  =  1  THEN  16270 
16280  ZH1126)  =  ZZ(3) 

16290  ZN(56)  =  INT  (ZNl  /  (ZH(55)  I  ZH(57)I):  SOTO  15300 
16300  HONE  :  VTAB  3:  PRINT  'SCALE  FACTORS:* 

16310  VTAB  6:  PRINT  '//  RESISTANCE  STANDARD  ';:ZZ(3>  *  ZH(51):ZZ(4>  *  .1:11(5)  *  1E9:ZZ(6>  =  0:  GOSUB 
10100:  IF  ZEKO)  *  1  THEN  16310 
16320  ZN(51>  =  11(5) 

16330  VTAB  10:  PRINT  '//  CAPACITANCE  STANDARD  *;:ZZ<3>  *  ZH(52):ZZ(4)  =  0:11(5)  =  1E3:ZZ(6)  *  0:  GOSUB 
10100:  IF  ZEKO)  =  1  THEN  16330 
16340  ZM1S2)  =  11(5) 

16350  VTAB  14:  PRINT  'IHPEDANCE  SCALE  FACTOR  ';:ZZ(3)  =  ZH(53):ZZ(4)  =  1-.ZZ45)  =  1E9:ZZ(6)  =  0:  GOSUB 
10100:  IF  ZEKO)  =  1  THEN  16350 
16360  ZH153)  =  ZZ(3) 

16490  60T0  15300 

16500  HONE  :  VTAB  3:  PRINT  'NEASUREHENT  NODE:' 

16510  VTAB  5:  CALL  -  958:  PRINT  'INPUT(S):  0=2/1*:  HTAB  14:  PRINT  '1=1':  KTAB  14:  PRINT  *2=2':  HTAB 
14:  PRINT  *3=(2-l)/l':  HTAB  14:  PRINT  '4=1/42-1)':  HTAB  14:  PRINT  '5=162' 

16520  VTAB  11:ZZ(3)  =  ZH4 127) : ZZ (4)  =  0:ZZ<5>  =  5:ZZ(6)  =  1:  GOSUB  10*90:  IF  ZEKO)  =  1  THEN  16520 
16525  ZM4127)  =  11(5) 

16530  VTAB  13:  PRINT  'OPERATOR:  0=UNITY':  HTAB  14:  PRINT  '1*JN':  HTAB  14:  PRINT  '2=1/JN':  HTAB  14:  PRINT 
*3=JN*2':  HTAB  14:  PRINT  '4=l/«*2' 

16540  VTAB  18:ZZ(3)  =  ZH(116):ZZ(4)  =  0:ZZ(5>  =  4:ZZ(6)  =  0:  GOSUB  10100:  IF  ZEKO)  =  1  THEN  16540 
16545  ZH4116)  =  ZZ(3) 

16550  VTAB  20:  PRINT  'COORDINATES  0«A/B':  HTAB  14:  PRINT  'I=R/THETA':  HTAB  14:  PRINT  '2=L06.R/THETA' 

16560  VTAB  25:11(5)  =  ZH(46):ZZ(4)  =  0:11(5)  =  2:ZZ(6)  =  0:  GOSUB  10100:  IF  ZEKO)  =  1  THEM  16560 
16565  ZH446)  =  ZZ(3) 

16570  60T0  15300 

16600  HONE  :  VTAB  3:  PRINT  'RETRIEVE  SET-UP  FILE* 

16620  GOSUB  16680: 11(9)  =  77:  GOSUB  11300:  GOTO  15000 
16640  HONE  :  VTAB  3:  PRINT  'SAVE  SET-UP  FILE  TO  DISK' 

16660  60SUB  16680:ZZ(8)  =  77:  GOSUB  11400:  GOTO  15000 

16680  X  =  FRE  (0):  GOSUB  11500:  PRINT  :  PRINT  'FILE  NAME  HILL  BE  PRECEDED  BY  DHSS(P)-’:  PRINT  :  INPUT 
•ENTER  FILE  NAME  *  DHSS(P)-'|ZDf(0):ZF«  *  'DHSS(P)-'  ♦  ZD440) : ZDf 40)  *  ZF»:  RETURN 
16700  HOHE  :  VTAB  3:  PRINT  'EXECUTE  SET-UP  FILE* 

16710  FOR  II  *  1  TO  6:  GOSUB  17000:  NEXT  II 


16720  II  *  14:  GOSUB  17000:  IF  ZH(14)  =  1  THEN  FOG  II  =  11  TO  13:  60SUB  17000:  NEXT  II 
16730  II  =  20:  GOSUB  16900 

16740  II  *  21:  GOSUB  17000:11  -  24:  IF  ZN<23>  >  0  THEN  II  =  23 

16750  GOSUB  17000:11  =  25:  60SUB  17000:11  =  26:  GOSUB  17000 

16760  FOG  II  *  27  TO  30:  FOG  JZN  =  1  TO  2:  GOSUB  17100:  NEXT  JZN:  NEXT  II 

16770  II  -  34:  60SUB  17000:11  *  35:  60SUB  17000:11  -  46:  60SUB  17000 

16780  ZH(45)  =  0201:  IF  ZHI127)  >  0  THEN  ZH(45)  =  0100:  IF  ZHI127)  >  1  THEN  ZHI45)  =  0200:  IF  ZH(127I  > 

2  THEN  ZN(45)  =  0100 

16790  11  -  45:  GOSUB  17000:1!  -  68:  60SUB  16900 

16800  JZN  =  0:11  =  85:ZH<85)  =  1:  60SUB  17200:  GEN  DISPLAY 

16802  JZN  =  1 : ZH (85)  =  0:  GOSUB  17200 

16G04  JZN  *  2: ZH(BS)  =  1:  GOSUB  17200 

16806  JZN  -  3:ZH(85)  -  1:  GOSUB  17200 

16810  JZN  *  1:11  «  86: ZH<86)  -  0:  60SUB  17200:  GEN  DISPLAY 

16811  JZN  -  2:ZN<86)  =  0:  GOSUB  17200 

16812  JZN  =  3: ZN(86)  =  ZHU16):  GOSUB  17200 

16813  JZN  =  4: ZH (86)  =  0:  GOSUB  17200 

16814  JZN  -  5: ZHI66)  *  0:  60SUB  17200 

16816  11  -  108:  60SUB  17000:11  =  109:  60SUB  17000 
16820  GOTO  1000 

16900  GESTOGE  :  FOG  I  -  1  TO  II:  GEAD  H4:  NEXT  I:H<  -  V  ♦  HI  ♦  GOSUB  21000:  GETUGN 

17000  GESTOGE  :  FOR  I  =  1  TO  11:  GEAD  HI:  NEXT  I:H* . .  NO  ♦  STGt  (ZN(ID)  ♦  GOSUB  21000:  GETUGN 

17100  FOG  JZN  =  1  TO  2:  GESTOGE  :  FOG  I  *  1  TO  II:  GEAD  Ht:  NEXT  I:H»  «  V  ♦  NO  ♦  STG*  (JZN)  ♦  *,*  ♦ 

STG4  (ZH(ID)  ♦  *’•;  GOSUB  21000:  NEXT  JZN:  GETUGN 

17200  GESTOGE  :  FOG  I  =  1  TO  II:  GEAD  NO:  NEXT  I:H4  »  *:’■  +  HO  ♦  STGO  (JZN)  ♦  *,*  ♦  STGO  (ZH(II))  ♦ 

*’*:  GOSUB  21000:  GETUGN 
18000  REN  CLOCK 

18010  ZZ (3)  *  PEEK  (  -  16384)  -  128:  GEN  KEY  PUSH 

18020  IF  ZZ  <3)  *  1  THEN  60T0  1000:  GEN  ABOGT 

18080  IF  ZZ (3)  «  19  THEN  18220:  GEN  STATUS 

18090  IF  ZZ(3)  =  20  THEN  ZH<47>  =  999:  GEN  TERHINATE 

18100  IF  ZZ (3)  =  21  THEN  ZH(47)  =  0 

18190  RETURN 

18200  REN  TINE 

18220  II  =  1:H*  =  •:,?THO,■:  GOSUB  24000:ZZ(3)  *  ZZ(II):H*  =  VTTHr*:  GOSUB  24000:ZZ(4)  *  ZZ(II):ZD$( 

0)  =  ’TINE  IS  1  ♦  STR4  (ZZ(3)I  ♦  V:H*  «  STG*  (ZZ(4)):  IF  ZZ( 4)  <  9  THEN  H*  =  ‘O’  ♦  H* 

18240  ZD4(0)  =  ZD*(0)  ♦  H$:  VTAB  1:  HTAB  28:  PRINT  ZD»(0):  RETURN  :  REN  READ  CLOCK 
19900  REN  CHAIN  SUBROUTINE 

19920  HONE  :  VTAB  9:NPI  =  1:  PRINT  ’LOADING  DHSS.*ZD*(0> 

19940  PRINT  CW*  <4);aBL0AD  CHAIN, A520,S*SBl;t,Dl,V*VB2 
19960  RETURN 
20000  ZZX  *  12 

20010  ZD*(4)  »  *8"  ♦  CHRI  (95)  ♦  CHR*  (ZZX  ♦  64):ZDf(5>  =  *8'  ♦  CHR*  (95)  +  ■?•  ♦  CHRt  (ZZX  ♦  32)  ♦ 

• 

20020  ZDK6)  *  CHR*  (4)  ♦  MN#0,:ZD»(7)  *  CHR*  (4)  ♦  *PR»0* 

20030  ZD* (8)  *  CHR*  (4)  ♦  *PR#3,:AA*  =  ,8*  ♦  CHR*  (95)  ♦  •  < ; ’ ■  s BBB  =  ■X'B’tfH:' 

20900  RETURN 

21000  GEN  SEND  DATA  TO  FRA 
21030  PRINT  ZD* (8):  PRINT  ZD*(5) 

21040  PRINT  N« 

21050  PRINT  ZD* (6) 


21470  PRINT  ZD$(7) 

21S00  RETURN 

24000  REN  SATA  RETRIEVE 

24030  PRINT  ZDt (8):  PRINT  ZD«(S) 

24040  PRINT  CHR4  (13):  PRINT  Ht:  60SUB  28500:  PRINT  ZDt<7):  PRINT  ZDt (6):  RETURN 

28500  PRINT  CHRI  (13);ZDI(4);':*:  PRINT  ZD*(7):  VTAB  23 

28501  ON  II  -  1  60T0  28504,28506,28508 

28502  INPUT  ZZ(1):  PRINT  ZD«<8):  RETURN 
28504  INPUT  ZZ(1),ZZ(2>:  PRINT  ZD! (8):  RETURN 
28506  INPUT  ZZ(1),ZZ(2),ZZ<3>:  PRINT  ZDt (8) :  RETURN 
28508  INPUT  ZZ(1),ZZ(2),ZZ(3),ZZ(4):  PRINT  ZDt (8):  RETURN 
29000  REN  KEITHLEY  HODEL  192  SET-UP 

29200  ZDt(O)  =  Nt  *  ■R0K0B0S3H0Z0N1 ' 

29220  VTAB  23 

29240  PRINT  ZDt<8):  PRINT  AAt;Z8t(0);BBt:  PRINT  ZDt(7):  INPUT  ZDt<0>:  PRINT  ZDt(6):  VTAB  23:  PRINT  *• 
29260  ZZ(3)  *  VAL  <  NIDI  (ZDt(0),5)):  RETURN 
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I 
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